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ABSTRACT 



Aims. In this paper we study the possibility of testing CPT symmetry with Cosmic Microwave Background (CMB) measurements. 
Methods. Working with an effective lagrangian of the photon with CPT violation X ~ p^AyF''^ which causes the polarization vectors 
of the propagating CMB photons rotated, we determine the rotation angle Aa using the BOOMERanG 2003 and the WMAP3 angular 
power spectra. 

Results. In this analysis we have included the newly released TC and GC (I < 450) information of WMAP3 and found Aa = -6.2+3.8 
deg at 68% confidence level. 

Conclusions. This result increases slightly the significance for the CPT violation obtained in our previous paper (Feng et al. 2006) 
Aa = -6.0 ± 4.0 deg (Itr). Furthermore we examine the constraint on the rotation angle from the simulated polarization data with 
Planck precision. Our results show that the future Planck measurement will be sensitive to Aa at the level of 0.057 deg and able to 
test the CPT symmetry with a higher precision. 
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1. Introduction 

In the standard model of particle physics CPT is a fundamental 
symmetry. Probing its violation is an important way to search 
for the new physics beyond the standard model. Up to now, CPT 
symmetry has passed a number of high precision experimental 
tests and no definite signal of its violation has been observed 
in the laboratory. So, the present CPT violating effects, if exist, 
should be very small to be amenable to the experimental limits. 

The CPT symmetry cou ld be d ynamically violated 
in the expanding universe dLi et al. 2 007). The cosmo- 
logical C PT violation mechanism investigate d in the ht- 
erature (iLi et al. 20021: L i & Zhang 2003; iLi et al. 2004t 



iFeng et al. 2005 : Li et al. 20071) has an interesting fea- 
ture that the CPT violating effects at present time 
are too small to be detected by the laboratory exper- 
iments but large enough in the early universe to ac- 
count for the generation of matter-antimatter asymmetr y 



(ILi et al. 2002t iLi & Zhang 2003t iLi et al. 2004t iLi et al. 2007h . 
And more importantly, this type of CPT violating effects could 
be accu mulated to be observable in the cosmological experi- 
ments (Feng et al 2001 iLi et al. 20071 iFeng et al. 20061) . With 
the accumulation of high quality observational data, especially 
those from the CMB experiments, cosmological observation 
becomes a powerful way to test the CPT symmetry. 

Here we study the CMB polarizations and CPT vi- 
olation in the photon sector with an effective lagrangian 
(ICarroU et al. 1990HCarroU & Field 199 iT) : 



(1) 



where Xcs ~ PfiAyF^'^ is a Chern-Simons term, is an ex- 
ternal vector and pi"'' = (l/2)e'"'^'^Fp^ is the dual of the 



electromagnetic tensor. This Lagrangian is not gauge invari- 
ant, but the action is gauge independent if dyPf^ - d^py. This 
may be possible if p/i is constant in spacetime or the gradi- 
ent of a scalar field in the quin t essential baryo-/lep togenesis 
dLi & Zhang 20031: iLi et al 2002t iDe Felice et al. 20031) or the 
gradient of a func tion of the Ricci scalar in gravita tional 
barvo-/leptogenesis dLi et al. 20041 : lOavoudiasl et al. 2004 ). The 
Chern-Simons term violates Lorentz and CPT symmetries when 
the background value of is nonzero. 

One of the physical consequences of the Chern-Simons term 
is the rotation of the polarization directio n of electromagneti c 
waves propagating over large distances dCarroll et al. 199()l) . 
From the Lagrangian ([1]), we can directly obtain the equation 
of motion for the electromagnetic field: 

V^(V^A'' - VA^) = -p^e^W(v^A, - V^Ap) . (2) 

After imposing Lorentz gauge condition V^A'' = 0, it becomes: 

VpVM" + RIA'' = -p^^-P-(T7pA^ - V^Ap) , (3) 

where RJ^ is the Ricci tensor. With the geometric optics approxi- 
mation, the solution to the equation of motion is expected to be: 
A'' - Re[(fl'^ + ebf + e^c^ + ...)e'^^'^], where e is a small num- 
ber. With this ansatz, one can easily see that the Lorentz gauge 
condition implies fc^a'' = 0, where the wave vector = Vp5 is 
orthogonal to the surfaces of constant phase and represents the 
direction which photons travel along with. The vector a'' is the 
product of a scalar amplitude A and a normalized polarization 

1 . Hence in the Lorentz gauge, 



vector e*^, a*^ 



: As^, with Sfjs'' 
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the wave vector kf^ is orthogonal to the polarization vector s''. 
Substituting this solution into the modified Maxwell equation (O 
and neglecting the Ricci tensor we have at the leading order of e 
the equation is kf^k^ = 0. It indicates that photons still propagate 
along the null geodesies. The effect of Chern-Simons term ap- 
pears at the next order, k^V^e'' - -pf,£^^'^°'kp£a-. We can see that 
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the Chern-Simons term makes A;'^V^,e'' not vanished. This means 
that the polarization vector e" is not parallel transported along 
the light-ray. It rotates as the photon propagates in spacetime. 

We consider here the spacetime described by spatially flat 
Friedmann-Robertson- Walker (FRW) metric. The null geodesies 
equation is {k^)^ - k'k' - 0. We assume that photons propagate 
along the positive direction of x axis, i.e. = A:' , 0, 0) and 
A;' = Gauge invariance guarantees that the polarization vector 
of the photon has only two independent components which are 
orthogonal to the propagating direction. So, we are only inter- 
ested in the changes of the components of the polarization vec- 
tor, and s^. Assuming = po to be a non-vanishing constant, 
we obtain the following equations: ds-^/dA + 'Hk^s^ - pok^^s^, 
ds^/dA + 'Hk°e^ = -pok^e^, where we have defined the affine 
parameter A which measures the distance along the light-ray, 
k'' = dx^/dA, and the reduced expansion rate "H = d/a. The 
polarization angle is defined as x = arctan(£''/e^). It is easy to 
find that the rotation angle is 



= XQ -Xz = - podri^ I 



Po 



dt 



where the subscript z is the redshift of the source when the light 
was emitted. For CMB photons, the source is the last scattering 
surface with z — llOOQ. The subscript indicates the present 
time. As we know, a vector rotated by an angle Ax in a fixed 
coordinates frame is equivalent to a fixed vector observed in a 
coordinates frame which is rotated by -A^. So, with the notion 
of coordinates frame rotation, the rotation angle is 



Aa = -Ax 



dt 

Pi) — = Por. 
a 



(5) 



with r, being the comoving distance of the light source away 
from us. This phenomena is known as "cosmological birefrin- 
gence". This rotation angle Aa can be obtained by observing 
polarized radiation from distant sources such as radio galaxies, 
quasars and CMB. 

The Stokes parameters Q and U of the CMB polarization 
can be decomposed into a gradient-like (G) and a curl-like (C) 
component ( Kamionkowski et al. 1997). For the standard theory 
of CMB, the TC and GC cross-correlation power spectra vanish. 
With the existence of cosmological birefringence, the polariza- 
tion vector of each photon is rotated by an angle Aa, and one 
would observe nonzero TC and GC correlations, even if they are 
zero at the last scatterin g surface. Denoting the rotated qu antities 
with a prime, one gets dFeng et al. 20051: [Lue et al. I999h : 



C;^"^ = C/'"sin(2AQ') , 
1 
2 

^'TC ^ cj^ cos(2Aa), 

' Besides the CMB photons which come from the last scattering sur- 
face, we might observed the different CMB photons which travelled 
different distances as well. However, we are only interested in linear 
perturbations of CMB photons in this paper. CMB polarizations are al- 
ready linear phenomena. They are not existent at the zeroth order. When 
calculating the variations of polarizations due to CPT violation in the 
perturbation theory up to linear order, we may ignore the fluctuations of 
the travelling distances of CMB photons. Otherwise, these fluctuations 
combined with polarizations would give higher order result which are 
beyond the scope of this paper. So, in this paper, we only consider linear 
perturbations and we can assume that each CMB photon detected by us 
travelled the same distance. 



-(Cf^ - Cf^) sm{4Aa) 
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Fig. 1. One dimensional distribution on the rotation angle Aa 
from WMAP3 and B03 data. 



(4) C, 



'GG 



C?^ cos^(2Aa') + C'f'- sin^(2AQ') 



cc , 



(j'cc ^ cf^ cos\2Aa) + Cf^ sm\2Aa) , (6) 
while the temperature power spectrum TT remains unchanged. 

2. Method and Results 

In our study we make a global analysis to the CMB data 
with the public available Markov Ch ain Monte Carlo pack- 
age CosmoMCQ ( iLewid & Bridle 2002 ^. which has been mod- 
ified to allow the rotation of the power spectra discussed above, 
with a new free parameter Aa. We assume the purely adia- 
batic initial conditions and impose the flatness condition moti- 
vated by inflation. Our most general parameter space is: P = 
(a>b, oJc, ®s, T, lis, log[10'"Ai], r, Aa), where lOh = Q.th^ and cjc = 
Q.ch^ are the physical baryon and cold dark matter densities rel- 
ative to the critical density, 0, is the ratio of the sound hori- 
zon to the angular diameter distance at decoupling, t is the opti- 
cal depth to re-ionization. A , and n , characterize the primordial 
scalar power spectrum, r is the tensor to scalar ratio of the pri- 
mordial spectrum. For the pivot of the primordial spectrum we 
set ksQ - O.OSMpc"'. In our calculation we have assumed that 
the cosmic rotation angle is not too large and imposed a conser- 
vative flat prior -7r/2 < Aa < 7r/2. 

In our calculation we combine the full d ata of WMAP3 
(Spergel et al. 2007; Page et al. 2007; Hins haw et al. 20071; 
iJarosik et al. 20071) (including the information of TC and GC 
power spectra (/ < 450)) and BOOMERanG 2003 (B03 ) 
(iJones et al. 20061; iPiacentini et al. 20061; iMontrov et al. 20061) . 
We calculate the likelihood of TT, TG, GG and CC power 
spectra using the routine for computing the likelihood supplied 
by the WMAP team. As for the TC and GC power spectra data 
(/ < 450), which are just the preliminary suboptimal results 
currently, we simply assume the Gaussian likelihood functiorQ: 



-Ctc/gc - exp 



^TC/GC 



" http://cosmologist.info/. 

' We are very grateful to Professor Eiichiro Komatsu (the member 
of WMAP group) for the email communications on the use of the TC 
and GC data of WMAP3. Because the current released TC and GC 
information of WMAP3 are not very accurate and do not include the 
correlations between different / values, our simple Gaussian likelihood 
function Eq.Q is accurate enough for the illustrative purpose. 
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Z 



^TC/GC _ ^TC/GC \2 



(7) 



where C^'xJ'^'^ and C^'^J'^'^ denote the theoretical value and the 

l,th l,obs 

experimental data of the TC or GC power spectra of WMAP3, 
o"! denotes the variance of estimated power spectra at each multi- 
pole. 

In the computation of the CMB spectra, we have consid- 
ered lensing contributions. Furthermore, we make use of the 
Hubble Space Telescope (HST) measurement of the Hubble pa- 
rameter Ho = lOOh km s ' Mpc ' by m ultiplying a Gaussia n 
likelihood function h = 0.72 + 0.08 CFreedman et al. 20011) . 
We also impose a weak Gaussian prior on the baryon density 
Q.hh^ ^ 0.022 + 0.002 (Icr) from the Big Bang Nucleosynthesis 
jBurles et al. 2001). Simultaneously we will also use a cosmic 
age tophat prior as 10 Gyr < fo < 20 Gyr. 

In Figure [T] we plot our one-dimensional constraint on the 
rotation angle Aa from the combined WMAP3 and B03 data 
and find that the current CMB polarization data favor a nonzero 
rotation angle of the photons. The best fit value of the rotation 
angle is Aa = -0.122 = -7.0 deg. Marginalizing the posterior 
distributions, we find that the mean value of the rotation angle 
is: 



Aa = -0.108 + 0.067 = -6.2 + 3.8 deg (Icr) , 



(8) 



which deviates the unrotated case Aa - more than Icr 
and gives a weak evidence for CPT violation. This re- 
sult is frequency independen t, while the Faraday Rotation 
jScannapieco & Ferreira 1997 1). which also can give nonzero TC 
and GC power spectra, is frequency dependent and of the much 
smaller order than ours. On the other hand, currently several high 
precision experimental tests have confirmed the CPT symme- 
try and do not detect apparent CPT violation in the laboratory. 
However, these two results are consistent. From equations ( |4|5t 
we can find that the size of the effects of CPT violation on CMB 
polarization power spectra is an integrated effect and enhanced 
during the propagation time, while the laboratory experiments 
just test the CPT symmetry at present. It means that the CPT 
violating effect should be too small to be detected at any specific 
time and enhanced during the propagation time to be detectable 
by the CMB measurements. 

In Figure |2] we illustrate the effects of the rotation angle Aa 
on the power spectra for two cases Aa - and Aa - -0.108. 
The basic cosmological parameters we choose for these plots 
are the best fit values of global fitting: Q/,/?^ = 0.0228, Q.,„h^ - 
0.1320, T = 0.099, //o = 73.1, n, = 0.965, A, = 2.2 x 10"^ 
and r - 0.00866. One can see from Figure 2 that the C-mode is 
very sensitive to the rotation angle Aa. Direct measurements of 
the TC and GC power spectra are crucial and would give more 
stringent constraints than ones of other power spectra. 

From the analysis of data points of TC and GC power spec- 
tra of WMAP3 and B03, we can find that this negative rotation 
angle is slightly preferred. Firstly, in the B03 data, the TC power 
at Z ~ 250 and / ~ 350 are both negative, while it is positive at 
I ~ 450. The GC power at I ~ 250, / ~ 350 and I ~ 450 are 
all negative. Based on the equations (|6]) and Figure |2l we can 
see that the TC and GC power spectra of B03 really help to ob- 
tain this negative rotation angle. On the other hand, the quality 
of TC and GC power spectra of WMAP3 seems not very good. 
The data points are distributed around the line Aa - propor- 
tional and the error bars are very large. Therefore, the WMAP3 
polarization data can not give much more contribution for test- 



TABLE I. Assumed experimental specifications. The noise parameters 
At- and Ap are given in units of /jK-arcmin. 



,/sky 


/max (GHz) 


^fwhm 


At- 


Ap 


0.65 


2500 100 


9.5' 


6.8 


10.9 




143 


7.1' 


6.0 


11.4 




217 


5.0' 


13.1 


26.7 



ing the CPT symmetry and not improve the previous result of 
rotation angle significantly. 

Recendy Feng et al. (2006) used the WMAP3 and B03 po- 
larization data to constrain the rotation angle and found the sim- 
ilar result: Aa - -6.0 + 4.0 deg at Icr confidence level. As we 
show above, the TC and GC information are very important for 
this type of analysis. In their analysis they did not include the 
TC and GC information of WMAP3 and set the TC and GC of 
power spectra of WMAP3 to be zero, Cj'^ - and Cf^ - 0. 
This factor should affect the final result. Comparing these two 
results, we can find that the mean value of our result moves to- 
ward direction away from the zero and the error bar shrinks a 
little, which strengthens the conclusion on the CPT violation. 

Moreover, Cabella et al. (2007) perform a wavelet analysis 
of the temperature and polarization maps of the CMB delivered 
by the WMAP3 experiment which include the information of 
TC and GC power spectra. They set a limit on the CMB photon 
rotation angle Aa = -2.5+3.0 deg (Icr). They found no evidence 
for the CPT violation from WMAP3 polarization data. However, 
in their analysis, they have not included the information of B03 
polarization data. We also made an analysis without the B03 data 
and found that Aa - -1.6 + 10.3 deg (Icr) which means that 
the WMAP3 data only did not give any significant evidence for 
CPT violation. Therefore, we believe that they should obtain the 
similar result with ours if they also include the B03 polarization 
data in their analysis. 

Obviously the current CMB polarization data are not good 
enough to verify this possible CPT violation. We need more ac- 
curate CMB data such as the Planck measurement^ in the near 
future. In Table I we list the assumed experimental specifica- 
tions of the future Planck measurement. The likelihood function 



is X oc e\p(-xlff/2) and 



(9) 



where /s^y denotes the observed fraction of the sky in the real 
experiments, A is defined as: 

A = Cf^iCf^Cf^ - (Cf^f) + Cj^iCf^Cf^ - Cj^Cf^) 
+ Cj'^iCj'^Cf'^ - Cj^C^^) + Cf "^(Cf "^cf^ - C]^Cf^) 



+ cf'^icj^cf'^ - (cf^f) + cf'={cj^c]'= - cYcf^) 

+ Cf^iCf^Cf^ - Cf^Cf^) + Cf^iCj'^Cj^ - Cj'^Cf^) 
+ Cf^iCj'^C'f'^ - (Cj^f) , 



(10) 



and |C| and \C\ denote the determinants of the theoretical and 
observed data covariance matrices respectively. 



|C| 



cj^cf^cf^ + ic]^c]^c^^ ■ 



■ cF(c?^f 



_rGC 



|C| = cj^cf'^cf'^ 



•cc, 

TCATCACC 



+ 2Cj ""Cj *- cy"- - c/ ' (C^^f 



-cf^icj'^f - Cf^iCj^f 



(11) 



http://sci.esa.int/science-e/www/area/index.cfm?fareaid=17/. 
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The likelihood has been normalized with respect to the maxi- 
mum likelihood - 0, where Cf^ - Cf^. For more details 
of this calc ulation we refer the readers to our previous compan- 
ion paper (IXia et al. 2007h . By the simulated data with Planck 
accuracy, we find that the standard deviation of the rotation an- 
gle will be significantly reduced to cr = 0.057 deg. This result is 
much more stringent than the current constraint. Assuming the 
mean value unchanged, the cosmological CPT violation will be 
confirmed around lOOcr confidence level with the future Planck 
measurement. 

In summary, probing the violation of fundamental symme- 
tries is an important way to search for the new physics beyond 
the standard model. In this paper we have determined the rota- 
tion polarization angle Aa with the CMB data from the WMAP3 
and BOOMERanG 2003. We find that a nonzero rotation angle 
of the photons is mildly favored: Aa - -6.2 + 3.8 deg (Icr). 
Our result shows a small violation of the CPT symmetry, albeit 
not conclusive with the present CMB data. With Planck data our 
simulations indicate that this type of CPT violation could be 
confirmed significantly, or CPT symmetry will be verified with 
a high precision. 
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